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Abstract 



The saturation of QCD chiral sum rules of the Weinberg-type is analyzed using ALEPH and OPAL 
experimental data on the difference between vector and axial- vector correlators (V-A). The sum rules 
exhibit poor saturation up to current energies below the tau-lepton mass. A remarkable improvement 
is achieved by introducing integral kernels that vanish at the upper limit of integration. The method is 
used to determine the value of the finite remainder of the (V-A) correlator, and its first derivative, at 
zero momentum: 11(0) = -4L M = 0.0257 ± 0.0003 , and ff'(0) = 0.065 ± 0.007 GeV -2 . The dimension 
d — 6 and d — 8 vacuum condensates in the Operator Product Expansion are also determined: 
< 6 >= -(0.004 ± 0.001) GeV 6 , and < 8 >= -(0.001 ± 0.006) GeV 8 . 
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Since the pioneering work of Shifman, Vainshtein and Zakharov [1], a few thousand papers have been 
published on applications of the QCD sum rule method in all corners of low energy hadronic physics. Un- 
avoidably results from different collaborations were not always consistent [2]. The main reason for these 
inconsistencies was frequently the impossibility of estimating reliably the errors in the method. With the 
advent of precise measurements of the vector (V) and axial-vector (A) spectral functions, obtained from 
tau-lepton decay [3]- [4], an opportunity was opened to check the precision of the QCD sum rules in the 
light-quark sector of QCD. In this note we would like to present a critical and conservative appraisal of 
chiral sum rules of the Weinberg type [5] , as they are confronted with experimental data for the spectral 
functions. This kind of sum rules involve the difference between the vector and the axial- vector corre- 
lators (V-A), which vanishes identically to all orders in perturbative QCD in the chiral limit. In fact, 
neglecting the light quark masses, the (V-A) two-point function vanishes like 1/q 6 in the space-like region, 
where the scale 0(300 MeV) is set by the quark and gluon condensates. In the time-like region the chiral 
spectral function pv-a(q 2 ) should also vanish for large Q 2 = —q 2 , but judging from the ALEPH data[3], 
the asymptotic regime of local duality may not have been reached in r-decay . Under less stringent as- 
sumptions one expects global duality to hold in the time like region; in particular, this should be the case 
for the Weinberg-type sum rules. Surprisingly, these sum rules also appear to be poorly convergent. A 
possible source of duality violation could be some non-perturbative contribution to the correlator (e.g. due 
to instantons) which falls off exponentially in the space-like region but oscillates in the time-like region. 
If the duality violations were due to this source, then there would be a simple recipe (introduced 30 years 
ago [6]) to improve convergence. 

In a previous publication [7] we studied some QCD chiral sum rules of the Weinberg type, and their 
saturation by the ALEPH data. In particular, we showed that a remarkable improvement of this satu- 
ration can be achieved by introducing a polynomial integration kernel which vanishes at the upper limit 
of integration. However, no detailed quantitative error analysis was performed in [7]. In this note we 
reexamine the saturation of several QCD chiral sum rules using the ALEPH [3], as well as the OPAL data 
[4], and paying particular attention to the error analysis. We obtain an updated determination of L w , 
the scale independent part of the coupling constant of the relevant operator in the 0(p 4 ) counter terms 
in the Lagrangian of chiral perturbation theory [8] . This quantity is related to the finite remainder of the 
(V-A) correlator at zero momentum. We also determine the finite remainder of the first derivative of the 
(V-A) correlator at zero momentum, which is related to the 0(p e ) counter terms. Finally, we introduce 
combinations of QCD chiral sum rules which allow for a determination of the (V-A) dimension d = 6 and 
d = 8 vacuum condensates. The former can be extracted with reasonable precision, while the latter is 
affected by much larger uncertainties. 
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We begin by denning the vector and axial-vector current correlators 

U^(q 2 )=i I d A xe^ <0\T(V,(x) Vj (0))\0 > 

= (-SV I 2 + %<1») n v(<? 2 ) , (1) 
n^(9 2 ) = i J d*x e** < 0|T(A„(aO 4(0))|0 > 

= (-sv i 2 + W") u a(i 2 ) - w» n (<? 2 ) , (2) 

where V^(x) =: q{x)^ )1 q{x) :, A fl (x) =: q(x)j ^ 759(0;) :, and q = (u,d). Here we shall concentrate on the 
chiral correlator Hy-A = IV — ^Ia- This correlator vanishes identically in the chiral limit (m q = 0), 
to all orders in QCD perturbation theory. Renormalon ambiguities are thus avoided. Non-perturbative 
contributions due to vacuum condensates contribute to this two-point function starting with dimension 
d = 6 and involving the four-quark condensate. The Operator Product Expansion (OPE) of the chiral 
correlator can be written as 



n(Q 2 )\v-A = £ c 2N+4 (Q 2 ^ 2 ) < o 2N+4 (^) > , 

N=l ^ 



(3) 



with Q 2 = —q 2 . It is valid away from the positive real axis for complex q 2 , and \q 2 \ large. Radiative 
corrections to the d = 6 contribution are known [9]. They depend on the regularization scheme, implying 
that the value of the condensate itself is a scheme dependent quantity. Explicitly, 



n(Q z )\v-A = 



32?r a s < qq > 2 



Q 6 



a s (Q 2 ) 

47T 



247 
~12 



0(l/Q 8 ) 



(4) 



in the anti-commuting 75 scheme, and assuming vacuum saturation of the four-quark condensate. Radiative 
corrections for d > 8 are not known. To facilitate comparison with current conventions in the literature 
it will be convenient to absorb the Wilson coefficients, including radiative corrections, into the operators, 
and rewrite Eq.(3) as 



00 _^ 



(5) 



N=l 



where we have dropped the subscript (V-A) for simplicity. We will be concerned with Finite Energy Sum 
Rules of the type 
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W(s ) = / ds f(s)p(s) , (6) 
Jo 

where f(s) is a weight function, and the hadronic spectral function p(s) = pv{ s ) — Pa(s), with pv,a{s) = 
^ImUv,A(s) (pion pole excluded from pa(s)). For instance, if /(s) = s N (N — 0, 1, 2, ...), then one obtains 



/ ds s N p(s) = fi5 N0 + (~) N < CW > (N = 0,1,2,...), (7) 
Jo 

where / w = 92.4 ± 0.26 MeV [10]. For N = 0, 1 Eq.(8) leads to the first two (Finite Energy) Weinberg 
sum rules, while for N = 2, 3 the sum rules project the d — 6, 8 vacuum condensates, respectively; notice 
that in the chiral limit < C 2 >=< O4 >= 0. To first order in a Sl radiative corrections to the vacuum 
condensates do not induce mixing of condensates of different dimension in a given FESR [11]. Wc shall 
also consider the chiral correlator, and its first derivative, at zero momentum; the finite remainder of these 
being given by the sum rules 



fi(0) = / - p(s) , (8) 

Jo s 

n'W - f° $ p{s) , (9) 

Jo s 

where p(s) does not contain the pion pole. Equation (8) is the Das-Mathur-Okubo (Finite Energy) sum 
rule [5] . The finite remainder 11(0) = — 4Lio , where L10 is a counter term of the 0(p 4 ) Lagrangian of 
chiral perturbation theory, can be expressed as 



11(0) = -4L 
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1 



-fl <rl> -F A 



0.026 ±0.001 



(10) 



where < r\ > is the electromagnetic mean squared radius of the pion, < r\ >= 0.439 ± 0.008 fm 2 [12], 
and Fa is the axial- vector coupling measured in radiative pion decay, Fa = 0.0058 ±0.0008 [10]. Similarly, 
II'(O) is related to the 0(p 6 ) counter terms. 



As mentioned earlier, the saturation of the various chiral sum rules can be considerably improved by 
introducing an integration kernel that vanishes at the upper limit of integration (s — so). We have tested 
a variety of such kernels searching for optimal saturation. The following results have been obtained using 
the ALEPH data for p(s), with the errors at each energy bin calculated from the error correlation matrix. 
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Use of the OPAL data [4] data leads to similar results, albeit with much lager error bands. Starting with 
the first Weinberg sum rule, Fig.l shows the left hand side of Eq.(7) for N = (curve(a)), together with 
the right hand side, i.e. / 2 (straight line (c)), as well as the modified sum rule (curve(b)) 

W 1 (s )= / ds (1 - -) p(s) . (11) 
Jo s o 

On account of the second Weinberg sum rule, curves (a) and (b) should be identical; the improved satu- 
ration achieved with Eq.(ll) being remarkable. Figure 1 can be used to present our criterion to judge the 
reliability of a QCD sum rule. The sum rule must be presented explicitly as a function of the upper inte- 
gration limit so- If the left hand side is a constant, then the spectral integral must also be approximately 
a constant, starting from 1 to 2 GeV 2 up to the maximum so of the data. From Fig.l we would extract 

fl = 0.008 ± 0.004 GeV 2 , (12) 

for curve (a), and 

fl = 0.0084 ± 0.0004 GeV 2 , (13) 

for curve (b), to be compared with the experimental value / 2 |bxp = 0.00854 ± 0.00005 GeV 2 . Curve (a) 
demonstrates the fact that if the spectral integral is not a constant then the experimental errors are quite 
irrelevant in a test of duality. It is very dangerous to pick up a small stability region to obtain a prediction 
(here one could choose the region around 2 GeV 2 ). 

In Fig. 2 we show Eq.(8) (curve(a)) together with the modified sum rule (curve (b)) 



fl + P* (1 _± 
so Jo s s ' 



11(0) = 2 ^ + - (1 - -Y p(a) . (14) 



From the optimized sum rule (14) we obtain the value (straight line (c)) 

n(0) = -4lio = 0.0257 ± 0.0003 , (15) 

which is considerably more accurate than the leading order chiral perturbation theory result, Eq. (10). 
The agreement between Eqs. (10) and (15) may be an indication that higher order chiral corrections to 
the Das-Mathur-Okubo sum rule are indeed very small. Figure 3 shows Eq.(9) (curve (a)) together with 
the optimized sum rule (curve (b)) 



n'(o) = ± n(o) - 3 % + / % (l - - 

s 4 Jo s s 



) 3 P(s) , (16) 
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the latter giving (curve(c)) 

fi'(O) = 0.065 ± 0.001 GcV~ 2 . (17) 

We turn now to the determination of the d = 6 and d = 8 vacuum condensates. In Fig. 4 we show < Oq > 
as obtained from Eq.(7) with N = 2 (curve (a)), together with the result from the improved sum rule 
(curve(b)) 

< O e >= -fl al + al ds (1 - -f p(s) , (18) 
Jo s o 

which gives 

< 6 >= -(0.004 ± 0.001) GcV 6 . (19) 
This result can be compared with the vacuum saturation expression 

< Oe > \ VS = -yJ a, | <qq>\ 2 ^ -1.1 x 10- 3 GeV 6 , (20) 

to leading order in a s , and where we used < qq >= —0.014 GeV 3 , and a s = 0.5, at a scale of 1 GeV. 
Radiative corrections increase this estimate by a factor of two. The result Eq. (19) confirms pioneer 
determinations from e+e~, as well as tau-lepton decay data [13]-[14] indicating that the vacuum saturation 
approximation underestimates the d = 6 condensate roughly by a factor of 2-3. 

Finally, for < C 8 > Fig. 5 (curve(a)) shows the result from Eq.(7) with N = 3, together with the improved 
determination from the sum rule (curve(b)) 

< O s >= 8 si fl - 3 st S(0) + si r - (1 - -) 3 (s + 3s Q ) p(s) , (21) 

Jo s s o 

which gives 

< O s >= -(0.001 ± 0.006) GcV 8 , (22) 

in the region where the condensate is approximately constant (sq ~ 1.75 — 2.5 GeV 2 ), and assuming, opti- 
mistically, that the stability region has been reached beyond sq ~ 2.5 GeV 2 . It should be clear from Fig. 
5 that no meaningful determination of < Og > is possible using the standard FESR, Eq.(7). As expected, 
with increasing dimensionality, i.c higher powers of s in the dispersive integrals, the accuracy of the deter- 
mination of the vacuum condensates deteriorates considerably. It should be noticed that the results (19) 
and (22) do not rely on the vacuum saturation approximation. They also include all radiative corrections, 
and are correct to first order in a s . At order o? s and beyond, there is no longer decoupling of condensates 
of different dimensionality in a given FESR [11]. However, one expects these higher order radiative cor- 
rections to the Wilson coefficients in the OPE to be small. There seems to be general agreement in the 
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literature on the size of the d — 6 condensate, but there exists a number of inconsistent QCD sum rule 
determinations of the value of the d = 8 condensate. The results range from < Og >= —(3.5 ± 2.0) x 10~ 3 
GeV 8 [15] to < 0$ >= (4.4 ± 1.2) x 1(T 3 GeV 8 [16]. Our result is consistent, within the larg e errors, with 
a recent determination [17] (this reference contains a detailed comparative study of the literature). 

The poor convergence of ordinary QCD chiral sum rules is rather intriguing, as one would have expected 
good saturation at relatively low energies, given the very rapid fall-off of the chiral (V-A) correlator (see 
Eq.(4)). However, extrapolating the chiral correlator from the space-like to the time- like region can pro- 
duce strong changes close to the real-axis. In fact, violations of local duality at the 100% level have been 
shown to be possible using realistic models of the heavy quark chiral correlator [18]. The remarkable 
improved saturation achieved by introducing weight functions that vanish on the real axis at s = so could 
be taken as an indication that although perturbative QCD works well in the space-like region, this may 
not be the case in the time-like region, or near the cut, at least at energies below so ~ 3.5 GeV 2 . Finally, 
by using the chiral (V-A) correlator we have been able to extract the value of the d = 6 vacuum conden- 
sate with reasonable accuracy; for the d = 8 condensate the result is affected by a large uncertainty. In 
contrast, were one to attempt a determination from the vector correlator, and separately from the axial- 
vector one, the results would be quite inconclusive. This is due to the very large current value of K-qcd 
(Aqcd — 400 MeV) which makes the perturbative QCD term in the OPE so big that it overwhelms the 
power corrections. We estimate that if Kqcd ^ 330 MeV then the FESR, and even Laplace transform sum 
rules, will be unable to provide a conclusive determination of the vacuum condensates. Earlier standard 
extractions of these condensates from electron-positron annihilation [13] and tau-lepton decay [14] relied 
on past values Aq C d — 100 — 200 MeV. With these values of K-qcd the perturbative QCD term in the 
OPE is dominant but not overwhelming, and the power corrections can be clearly discerned. 

We would like to conclude with a general comment. Mathematically, the extraction of QCD parameters 
from experiment via sum rules constitutes a so called ill posed inverse problem (analytic continuation of 
an imprecisely known function). Small changes in the input data lead to large changes in the output. The 
problem is stabilized by extracting only a small number of parameters. Given the present accuracy of the 
r-decay data, we conclude from our analysis that only the condensate < O e > can be extracted with some 
degree of confidence, and only a rough idea of the order of magnitude of < O g >can be obtained. This sit- 
uation cannot be remedied by mathematical tricks like employing Laplace or Gaussian integration kernels. 
Only with forthcoming more accurate data, can one expect to extract higher dimensional condensates. 
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Figure Captions 



Figure 1. Curve (a) is the standard first Weinberg sum rule, Eq.(7) with TV = 0, curve (b) is the modified 
sum rule Eq. (11), and curve (c) is the experimental value of 

Figure 2. The chiral correlator at zero momentum, 5(0), from the standard sum rule Eq. (8) (curve (a)), 
and from the modified sum rule Eq.(14) (curve (b)), the latter leading to the prediction Eq. (15) (curve(c)). 

Figure 3. The first derivative of the chiral correlator at zero momentum, fl'(0), from the standard sum rule 
Eq.(9) (curve (a)), and from the modified sum rule Eq.(16), the latter leading to the prediction Eq.(17) 
(curve (c)). 

Figure 4. The dimension-six vacuum condensate from the standard sum rule, Eq. (7) with N — 2 (curve 
(a)), and from the modified sum rule Eq.(18) (curve (b)). 

Figure 5. The dimension-eight vacuum condensate from the standard sum rule, Eq. (7) with N = 3 (curve 
(a)), and from the modified sum rule Eq.(21) (curve (b)). 
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